Study of W-exchange Mode 




Dong-Sheng Du, Ying Li^, Cai-Dian Lii 
Institute of High Energy Physics, P.O.Box 918(4), Beijing 100049, China 

We calculate the branching ratio of rare decay (f>I^ using the perturbative QCD factoriza- 

tion approach based on /ct factorization. Our result shows this branching ratio is (8.7 ± 1.4) x 10"'', 
which is consistent with experimental data. We hope the CLEO-C and BES-III can measure it more 
accurately, which will help us understand QCD dynamics and D meson weak decays. 
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The precise estimation of the branching ratio for the hadronic D decays is very important both in theoreti- 
cal side and experimental side. As a rare decay, — > (f)K'^ plays an important role in testing QCD dynamics 
and searching for new physics. In the standard model picture, this decay mode is a pure annihilation type 
decay, also called W^-exchange mode. It has been discussed with a large branching ratio to explain the big 



difference of lifetime between D'^ and D+ 



In the 80 's of last century, this decay's branching ratio has 



Br(D° ^ = (9.4± 1.1) X 10"^ (1) 

To our knowledge, this decay mode is never calculated quantitatively in QCD. So, it is necessary for us to 

reanalyze this decay mode seriously. 

Based on the factorization hypothesis, many D decays have been calculated in naive factorization approach 

II I 

and QCD factorization approach 5]. However, Annihilation diagrams can not be calculated easily for its 
endpoint singularity. It is well known that perturbative QCD (PQCD) factorization approach is successful 
in calculating two-body B meson decays The end point singularity can be regulated by Sudakov form 
factor and threshold resummation through introducing the transverse momentum fcy of valence quarks. 
Thus, PQCD approach can give converging results and have predictive power. Using this approach, people 
have calculated many B meson pure annihilation type decays quantificationally, and the results agree with 
data well 7 1 . 
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In this paper, we will calculate the (pK^ decay in PQCD approach. In this decay, the W boson 

exchange induces the four quark operator sc ud, and the ss quarks included in cf), are produced from 
a gluon. This gluon can attach to any one of the quarks participating in the four-quark operator. In the 
rest frame of D meson, the produced s and s quark included in final states have momenta of 0{Pk/2) and 
0{P^/2), respectively. Therefore the gluon producing them has momentum q ~ 0{Md/2), which is nearly 
1 GeV. This hard gluon can be treated perturbatively. Therefore the hard part calculation involves not only 
the four quark operator but also the hard gluon connecting ss quark pair. The factorization here means the 
six quark hard part calculation factorize from the non-perturbative hadronization characterized by meson 
light cone wave functions. 

We work at the D meson rest frame for simplicity. In light-cone coordinate system, the momentum of the 
D, (j) and meson can be written as: 

Pi=Pz5 = ^ (l,l,Oi),P2=P0 = ^ (l,r2,0^),P3=P;^o = ^ (0,l-r2,0x), (2) 

where r = M^jMj^ and we neglect the meson mass -M-^o compared with the large D meson mass. 
Because this decay is Z) — > PV mode, the transverse polarization of <^ meson gives no contribution. The 
longitudinal polarization vector of 4> meson is given as: 

Denoting the light (anti-)quark momenta in £), </> and mesons as fci, ki, and ^3, respectively, we can 
choose 

fci = (a;iP+,0,kiT), fc2 = (a;2P2+,0,k2T),fc3 = (0, X3P3-, k3T). (4) 
Then, integrating over , ^ and k^ , we get the decay amplitude: 

Amplitude ~ ^dxYdxidx-^hYdhxhidh-Jb^dhz 

Tr[(7(t)$B(a;i, 6i)$^(a;2, 62, e)$if (a;3, hi)H{xi, e, t)St{xi) 6"^^] , (5) 

where bi is the conjugate space coordinate of kj^, and t is the largest energy scale in ff, as the function 
in terms of Xi and bi. The large logarithms (Inmw/t) coming from QCD radiative corrections to four 
quark operators are included in the Wilson coefficients C{t). The large double logarithms (In^ajj) on the 
longitudinal direction are summed by the threshold resummation, and they lead to a jet function St{xi) 
which smears the end-point singularities on Xi. The last term, e~^^*\ contains two kinds of logarithms. One 



of the large logarithms is due to the renormalization of ultra-violet divergence In tb, the other is resummation 
of double logarithm In'^ b from the overlap of coUinear and soft gluon corrections. This Sudakov form factor 
suppresses the soft dynamics effectively. Thus it makes perturbative calculation of the hard part H reliable, 
is the wave function which describes the inner information of meson M. 
For this decay, the relevant weak Hamiltonian is 1^: 



(6) 



where Ci{^){i — 1,2) is the QCD corrected Wilson coefficient at the renormalization scale and the four 
quark operators Oi and O2 are 



Oi = {Sidi)v-A{UjCj)v-A,02 = {SiCi)v-A{Ujdj)v-A 



The four lowest order Feynman diagrams of (pK^ in PQCD approach are drawn in FIG Q according 

to this effective Hamiltonian. 





FIG. 1: Leading order Feynman diagrams of — * (fiK 

By calculating the hard part H at the first order of as, we get the following analytic formulas. With the 
meson wave functions, the amplitude for the factorizable annihilation diagrams in Fig^a) and (b) results 
in Fn as: 



^0 "'0 



Ef{tl)hi{x2,X3,b2,b3){{x3 - l)(l)4,ix2)(t)K{x3) 

2r{x3 - 2)rK(t)''^{x2)(f>K{x3) - 2x3rrif(/)^(a;2)(/)^(a;3)} - {{x2 - 1)00 (0:2 )(/i^(a;3) 

+ 2x2rrK4>^{x2)4'Ki^3) - '^r{x2 - 2)rif00(x2)0x(a;3)}-B/(i^)/i2(a;2, a;3, 62, 63) 



(8) 



In the function, Cp ~ 4/3 is the group factor of SU(3)c, and vk ~ mQK /Md, with took — toJ^/to^ (tos is 



the s quark current mass). The functions Ef is : 

Ef{t) = + C2(0/3)a,(i) e-^*(*)-^-W. (9) 

The hard scale Vs in the amphtudes are taken as the largest energy scale in the hard part H in order to kill 
the large logarithmic radiative corrections: 

1 1 ^ 9 ... n 1 1 



tl = T^la.yi{MD^Jl-x^.+X3.r^, — , —),tl = max(Mi5 VI - a;2, — , — ). (10) 

The functions hi{i = 1, 2) in the decay amplitudes consist of two parts: one is the jet function St{xi) derived 
by the threshold resummation, the other is the Fourier transformation of the propagator of virtual quark 
and giuon. They are given as: 
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hiix2,X3,b2,b3)^Stil-X3) i^-j HI,'\Md^{1 - X2){1 - X3 + x^r^) &2) 

X \^H^^\Md^JI - X3 + x^r^ b2)MMD^Jl - + X3r2 63)0(62 - 63) + (&2 ^ 63)} , (11) 

h2{x2,XiMM)^ St{l-X2) (^y^ iJ^''(A/l5 V(l - •^2)(1 - X3 + x^r^) 62) 
X {i/<'^(A/z3VT^62)Jo(Afo\/r^ 63)^^(62 - 63) + (62 ^ 63)} . (12) 

The amplitude for the nonfactorizable annihilation diagrams in Fig^c) and (d) results in 

Mg — 6AttCfM^ / dxidx2dx3 / 6id6i 62(^62 0D(a;i, 61) 

X {(2^3 - 1) (l}4>i^2)<t>Ki^3) - r {x2 ^X'i- 2)rK(l)l{x2)(f>K{x3) -r{x2- X3) rif0^(x2)0^(a;3) 
- r (a;2 + X3 - 2) rK(t)li,{x2)(l)K{x3) - r {x2 - X3) r^f (/>^(x2)0]^(a;3)}i;™(t^)/i*,^)(a;i, a;2, a;3, 61, 62) 
- {{X2- 1)(I)4,{x2)cI)k(x3) - r {x2 + X3 - 4) r;<- 0^(2:2) 0]^(x3) - r {-X2 + X3) rK(l)l{x2)(f>K{x3) 
-r{x2- X3) rK(l)%{x2)4>]i{x3) + r{x2+ X3)rK(l)%{x2)(j)]i{x3)]E,n{tl^)h^a\xi,X2,X^,bi,b2) , (13) 

where Xi dependence in the numerators of the hard part are neglected by the assumption xi <ti X2,X3. In 
the above Equation, some functions are defined as: 

E„,{t) = C2(i)as(t)e-^"(*)-^*(')-^'^('); (14) 
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ti, = max(Mz3y'|Fg.)|, Afo V(l ^ a^2)(l ~ X3 + X3r^), -, -); (15) 



hl;'Hxi,X2,X3,bi,b2) = 



I 



Ko{MDF^j)bi), for i^f,)>0' 



[fR^^\Mn,J\F^\b,), for ^^2, 



(16) 



with: 



Fii) = + 2^2 + 1)(1 - 2:3 + X3r2), ) = a;3(a;2 - Xi)(l - r^) - 1. (17) 

The total decay amphtude for (f)K'^ decay is given a.s A = fuFa + Ma- The decay width is then 

r(Z?o cflf) = ^^(1 - r^)\V:,V.sA\'. (18) 




FIG. 2: The D meson distribution amplitude 



For numerical analysis, we use the following input parameters 



0: 



fo = 230 MeV, fx = 160 MeV, = 241 MeV, mox = 1-7 GeV, 
\Vud\ - 0.9734, 114.1 = 0.974, T^o = 4.1 x lO-^^s^ /J ^ 220MeV. 



(19) 



The branching ratio obtained from the analytic formulas may be sensitive to many parameters especially 

those in the meson wave functions. The light (j) and K mesons' distribution amplitudes up to twist-3 have 

. . 1 . 

been used for many times in B decays 6|, and we don't list them here. The meson (0, K) wave functions 

are well constrained by hadronic B decays experiments. 

However, the heavy meson wave functions are still in discussion especially for D meson. In this wotIc, the 
D meson distribution amplitude we used has only one parameter ujjj, which is similar to B meson 



(j)Dix,b) = Nx'^{l - a;)^exp 



2uj' 



2- 2^^Dby 



(20) 



where iV is a normalization factor. For B meson, the peak appears at x — 0.1, because b quark is much 
heavier than the hght d quark. For D meson, the ratio of heavy (c) and hght quark (u) mass is rather 
smaller than that of B meson. So we adjust the parameter cj/^ = 0.5 GeV, which makes the distribution 
amplitude peak at x = 0.3. The shape of the D meson distribution amplitude is shown in FIG. 2. Using 
this wave function, we can get the form factor of ^ if as F^^^ = 0.80. This result is consistent with 
the previous calculation F^^'^ = 0.78 ± 0.04 . In fact, the heavy meson distribution amplitude can be 
eventually determined by radiative leptonic D meson decavsjioj. 

Here if we let ujd vary from 0.45 to 0.55, the branching ratios of D° (pK^ decay is: 



Br(D° ^ 0if") = (8.7±1.6) X 10"^ (21) 

which is consistent with the experimental measurement shown in eq.(^. From this calculation, wc find that 
the branching ratio becomes large when tj/j arise. Many other parameters such as Chiral breaking scale rriQK, 
CKM matrix elements also have large uncertainties, which will also enlarge the theoretical uncertainties for 
branching ratios. 

In general, if a process happens in an energy scale where there are many resonance states, this process 
must be seriously affected by these resonances. This is a highly non-perturbative strong interaction effect. 
Near the scale of D meson mass many resonance states exist, which may give large pollution to D decays 
calculation. So the final states interaction (FSI) may be important. However, we cannot calculate the FSI's 
contribution from the first principle. Although many people have discussed final states interaction in D 
decay ^|, there is large uncertainty in the calculation because they are usually model dependent. In our 
PQCD calculation we factorize the non-perturbative effects in meson wave functions, but neglect the soft FSI 
effect. Our numerical result is consistent with the experimental data well. It is a hint, that the contribution 
from soft FSI may not play an important role in this decay. We think the CLEO-C and EES can measure 
this VF-change channel more accurately, which will afford help for us to understand the dynamics of D 
decay. And the results also help us determine the D meson distribution amplitude. Since there is only tree 
operators contributing to this decay (only one kind of weak phase), there is no direct CP violation in the 
standard model. Any non-zero measurement of direct CP in this decay will be a signal of new physics. 

In a summary, we calculate the branching ratio of 13° (f)K^ in the perturbative QCD factorization 
approach without considering final states interaction. Our result indicates this branching ratio is very large 
comparing with other pure annihilation decay channels, because there is no Cabibbo suppression. This 



7 

branching ratio is about 8.7 x 10~^, and has been measured by CLEO. We hope the CLEO-C and BES-III 
can measure it more accurately, which will help us test QCD dynamics. 
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